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Abstract 

The deviation of primordial Helium production due to a variation on the 
difference between the rest masses of the nucleons is presented. It is found an 
upper bound 5(M n — M p ) < 0.129 MeV, between the present and nucleosyn- 
thesis epochs. This bound is used to analyze Wesson's theory of gravitation; 
as a result, it is ruled out by observation. 
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The so-called Hot Big Bang model provides a consistent description of the evolution of 
the Universe. This model depends on a set of universal parameters known as fundamental 
constants. However, a great number of theories of gravitation have appeared in which some 
of these constants do vary with time. In particular, scalar-tensor theories of gravity [1|-|3| 
make predictions which are in complete agreement with present experimental data while 
predicting a variation of Newton's gravitational constant over cosmological time scales. The 
physical embodiment of these theories allows a natural generalization of General Relativity 
and thus provides a convenient set of representations for the observational limits on possible 
deviations of Einstein's theory, making them a profitable arena for cosmology. 

The hypothesis that the gravitational interaction has changed over the history of the 
Universe (i.e. the gravitational parameter G or the rest masses of elementary particles 
depend on the time t) can be also analyzed in the framework of a 5-dimensional cosmology, 
proposed by Wesson 0. This theory is founded either on dimensional analysis |5],|6| as well 
as on reinterpretation of the five dimensional vacuum equations JF||| . A consequence of this 
relation is that the rest mass of a given body varies from point to point in space-time, in 
agreement with the ideas of Mach |5,IU|. This is a definite and testable prediction, especially 



when time intervals of cosmological order are considered P|,|TT 

In this letter, we find an upper bound to the variation of the masses of the nucleons over 
cosmological time intervals, from a comparison of the observed primordial abundance of 4 He 
with the theoretical variation induced by a changing mass. Subsequently, this upper bound 
is compared with the prediction of the 5-dimensional Wesson's theory of gravitation. 

If we consider that a variation of the rest masses of the particles had occurred between 
the epoch of primordial nucleosynthesis and ours, we can compute the deviation in the 4 He 
production from the Hot Big Bang model prediction due to this fact. The method we are 



going to apply is a generalization of the calculation made by Casas, Bellido and Quiros [12 
to fix nucleosynthesis bounds on the variation of the gravitational constant in Jordan-Brans- 
Dicke theory of gravity. The same method was recently used to limit other scalar-tensor 



theories with more general couplings functions [pl|. We know that - in the Hot Big Bang 
model - the primordial 4 He production is given by 

where A = exp(— (t nuc — £f)/r) stands for the fraction of neutrons which decayed into protons 
between t f and t nuc , with t { (t nuc ) the time of freeze out of the weak interactions (nucleosyn- 
thesis) |l4j| ) r the neutron mean lifetime, and x = exp(— (M n — M p )/kT) the neutron to 
proton ratio. If we make a variation on the rest masses of the particles, i.e. if we consider 
that a difference between the rest masses in the present and nucleosythesis epochs does 
exist, it is straightforward to compute an expression for the deviation in the 4 He primordial 
production, that reads 



$Y p = Y p \n [^-1 



AQ {Z) 



where we have defined AQ = M n — M p . 

We must also take into account that a variation in the rest masses of the particles will 
affect the neutron lifetime. The latter fact was not considered above since (0) represents 
only the explicit derivative with respect to AQ. A calculation of how a variation of the 



neutron lifetime affects the prediction on primordial 4 He has been already done in fl2,15 
It is given by, 

5Y P = 0.185— (3) 

T 

Noting that the dependence of r upon the masses is r oc Gp 2 AQ~ 5 oc AQ~ 5 , where 
is the Fermi constant and Mw is the mass of the bosonic mediator of weak interactions, it 
is easy to obtain, 

5t 5AQ 

T = AQ < 4) 

Thus, any variable rest mass theory will predict a primordial Helium abundance given 

by, 



Y 



p,var— rest— mass 



^ P ,std + SY P 



(5) 



where Yp jS td is the value predicted by the standard big-bang nucleosynthesis theory. 

We can find an upper bound for SY p summing up the two deviations referred above, i.e. 
equations (PD and (||), and comparing with the observed value Y t, s , 



SYJ < \Y ohs -Y std \+e<a 



(6) 



where e is an estimate of the observational error and a includes also estimates of theoretical 

errors. From [|l|,|ll| we estimate a < 0.01. But, it may also be possible that, due to small 

changes in nucleon masses, small changes in nuclear cross sections do occur. Since the 

functional dependence of cross sections with masses is generally unknown, we shall take into 

account these changes by arbitrarily doubling the theoretical error. Thus, we obtain 

5(AQ) 



AQ 



< 10% 



with AQ ~ 1.294 MeV 0, or equivalently 

5(AQ) < 0.129 MeV 



(7) 



(8) 



At this stage, we must work out the cosmological solution for the radiation dominated 
era. We shall consider the warped product M 4 x R, where M 4 is the ordinary 4-dimensional 
spacetime manifold and R correspond to the extra dimension. This leads to the line element 

A\t) 



ds^ = -dr + 



dx\ + dx\ + dx\ ) + e c{t) dx 



(9) 



(1 + A;(x2 + x2 + x2)/4) 2 
where A 2 {t) is the expansion scale factor and must be associated with the mass scale 
factor. This could be simplified considering a flat spatial section in M 4 (i.e. k = 0). This 
assumption is justified when one compares the order of magnitude of the different terms in 
the Einstein equation evaluated in the radiation era [19 



The suitable generalization of Einstein equation for the theory can be written as[| 



1 In what follows, latin indices A,B,..., run from to 4, greek indices from to 3 and latin indices 
from l to 3. 
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Gab = 8 it G T A b 



(10a) 



T AB = diag(p, -p, -p, -p, 0) 



(10b) 



with p (p) the density (pressure) of the radiation fluid, and G is Newton's gravitational 
constant. The equation of state is p = l/3p. The boundary condition that ought to be 
imposed is a smoothly matching at t = t eq ( equivalence time) with the dust-filled solution 
obtained in a previous paper |20|| . It is important to stress that in our previous work we 
required that both ( and d(/dt vanish at t = today. With these conditions, the masses 



of the fundamental particles can be set to their present experimental value [[T(J, and mass 



variations are negligible in short time scales (see eq.(|Tj)), which is consistent with the bound 
\fh / m\today ~ 10~ 12 yr^ 1 [^TJ. However, our results are to a large extent independent of the 
epoch in which the initial conditions for ( and d(/dt were imposed. 
The radiation dominated era solution is given by 



A 2 (t) = 2(3t 



(11a) 



C(t) = 2 In 



-1/2 



+ 




(lib) 



P{t) 



3 3 
+ 
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to 



(lie) 



(with (3 a constant)^. 



2 This can be seen by comparing the prediction of this theory for rh/m = Q/2 at t = today in the 
case where the initial conditions for £ and its derivative were imposed for instance at the Earth 



formation epoch, with the current experimental limits on rh/m [21|. 
3 This solution was previously obtained by Mann and Vincent [|l7j imposing different boundary 



conditions. It was also obtained by Gr0n [18|; from his work it is clear that the rate of change of 



the fifth metric coefficient depends on initial conditions. 



The way in which the mass is introduced in this formalism has an inherent ambiguity. 
The proposal of Ma 

"*(*) = 7? y/g^dx 5 (12) 

or, in the case of an xs-independent metric, 

c 2 

m(r) = pjv^F AZ (13) 

(A/ is the -finite- "length" of the body in the x$ direction 0), does not specify the tensorial 
character of the mass, which is implicitly introduced into the theory. In other words, there 
is no conclusive reason to mantain the covariant form of under the square root, and so 
m(t) can be also scale as 

m(r) = ^^Al (14) 

The bottom line in this idea is that the dimensional analysis used to define the relationship 
between the mass and the extra dimension x 5 = c 2 m / G does not exclude the covariant 
formulation £5 = c 2 m / G. We shall adopt hereafter the most favorable definition for the 
theory. 

Hence, we can estimate the variation of the mass from eq. (|llb[) . Since to ^> t cq ^> t nuc 
(to is the present age of the universe), we find 



C„ uc ^21n^/-^^21nl0 8 (15) 



Defining the quotient Am/m as 



Am = m(t n J-m(t ) = e _ c _ /2 
rriQ m(t ) 



we are able to see that the theory predicts 

SAQ 



100% (17) 



AQ 

which is in disagreement with the previous bound (0). Using a covariant formulation for the 
mass scale factor in ([T6| ) one would get an even worse disagreement. 
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Thereupon, the assumed relation connecting the 5 dimension and the particles rest 
masses is false, at least in the case we have explored. More general (i.e. x 5 dependant) 
metrics deserve more thorough analysis. We hope to report on this issues in a forthcoming 
work. 

We wish to thank Santiago Perez Bergliaffa and Harald Soleng for the critical reading 
of the manuscript and valuable comments. Remarks by an anonymous referee are gratefully 
acknowledged. This work has been partially supported by CONICET, UNLP and Fundacion 
Antorchas. 
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